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Abstract 

Positive linear systems on arbitrary time scales are studied. The the- 
ory developed in the paper unifies and extends concepts and results known 
for continuous-time and discrete-time systems. A necessary and sufficient 
condition for a linear system on a time scale to be positive is presented. 
Properties of positive reachable sets are investigated and characterizations 
of various controllability properties are presented. A modified Gram ma- 
trix of the system is used to state necessary and sufficient condition of 
positive reachability of a positive system on an arbitrary time scale. 

1 Introduction 

In many applications the variables that appear in a mathematical description 
take only positive or nonnegative values. Examples of such systems can be found 
in and [TT|, where also a theory of linear positive systems was developed. 
Usually the systems that are studied fall into two separate classes: continuous- 
time systems and discrete-time systems. In [llj all the problems are studied 
twice in these two separate settings. The characterizations of properties of 
positive systems for these two classes are sometimes similar, or even identical, 
and sometimes essentially distinct. 

Stefan Hilger in his Ph.D. thesis [TU] started the most successful attempt 
to unify the theories of continuous-time systems and discrete-time systems into 
one theory. It is based on the concept of time scale and the calculus on time 
scales. A time scale is a model of time. Time may be discrete or continuous, 
or partly continuous and partly discrete. The concepts of standard derivative 
used in the case of continuous time and forward difference used in the discrete 
time are unified into one concept of delta derivative. This allows to consider 
delta differential equations on arbitrary time scales. They generalize standard 
differential equations and difference equations. Theory of dynamical systems 
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on time scales was developed in [3]. Special attention was paid to linear delta 
differential equations. 

The interest in control systems on time scales dates back to 2004. The first 
results have concerned controllability, observability and realizations of linear 
constant-coefficient and varying-coefficient control systems with outputs [I] , [2] • 
Since then the literature on control systems on time scales has been rapidly 
growing, including also nonlinear systems. 

Controllability of continuous-time and discrete-time linear positive systems 
has been a subject of research since late 1980's [7], [E], [13]) [HI- Discrete- 
time systems appeared to be easier to deal with and it seems that positive 
cotrollability of such systems is now fully understood (see e.g. g], [5], [T2"]). 
On the other hand, only recently it was discovered that positive reachability 
of continuous-time systems requires very restrictive conditions to be satisfied 
[B], [TS]. Thus criteria of positive reachability for discrete-time systems and 
continuous-time systems are essentially different. 

In this paper we study linear positive constant-coefficients systems on arbi- 
trary time scales. The results presented here unify and extend corresponding 
results obtained for linear positive continuous-time and discrete-time systems. 
We prove necessary and sufficient conditions for a linear system x A = Ax + Bu 
on a time scale T to be positive. They involve the matrices A and B and the 
graininess function of the time scale, which describes distribution of the in- 
stances of time. We also study two controllability properties of positive linear 
systems: positive accessibility and positive reachability. Accessibility appears 
to be a property whose characterization does not depend on time scale. It is 
equivalent to standard controllability and expressed with the aid of the Kalman 
controllability matrix. As the criteria for positive reachability are completely 
different for continuous-time and discrete-time systems we have tried to develop 
methods which would result in the same statements for different time scales. 
We introduce a modified Gram matrix for a system on time scales, for which we 
select columns of B and choose different sets of integration for different columns. 
We prove that the system is positively reachable on an interval if and only if 
such a Gram matrix is monomial, i.e. each its column and each its row contain 
exactly one positive element. Then we show that from this characterization 
we can deduce known criteria for positive reachability of continuous-time and 
discrete-time systems. We also show that on nonhomogeneous time scales many 
properties known to hold for homogeneous time scales are no longer true. 

In Section [2] we recall basic material on positive systems, time scales and 
linear systems on time scales. Section [3J is devoted to positive control systems 
and Section [4] to positive accessibility and positive reachability. 

2 Preliminaries 

We introduce here the main concepts, recall definitions and facts, and set no- 
tation. For more information on positive continuous-time and discrete-time 
systems, the reader is referred to e.g. [9], and for information on time scales 
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calculus, to e.g. [3]- 

2.1 Positive math 

By R we shall denote the set of all real numbers, by Z the set of integers, and 
by N the set of natural numbers (without 0). We shall also need the set of 
nonnegative real numbers, denoted by R+ and the set of nonnegative integers 
Z+, i.e. N U {0}. Similarly, R h + will mean the set of all column vectors in 
R fc with nonnegative components and R* xp will consist of k x p real matrices 
with nonnegative elements. If A G R^ xp we write A > and say that A is 
nonnegative. A nonnegative matrix A will be called positive if at least one of 
its elements is greater than 0. Then we shall write A > 0. 

A positive column or row vector is called monomial if one of its components 
is positive and all the other are zero. A monomial column in R" has the form 
aek for some a > and 1 < k < n, where e& denotes the column with 1 at 
the kth position and other elements equal 0. Then we say that the column is 
k-monomial. An n X n matrix A is called monomial if all columns and rows of 
A are monomial. Then A is invertible and its inverse is also positive. Moreover, 
we have the following important fact. 

Proposition 2.1. A positive matrix A has a positive inverse if and only if A 
is monomial. 

It will be convenient to extend the set of all real numbers adding one element. 
It will be denoted by oo and will mean the positive infinity. We set R := RU{oo} 
and R + := R + U {oo}. If a G R then we define a + oo = oo. Moreover, for a G R 
and a > we set a/0 = oo and a/oo = 0. Of course oo > 0. If a matrix A 
has elements from R, then the notions of nonnegativity and positivity have the 
same meanings as before and are denoted in the same way. Addition of such 
matrices is defined in the standard way, but we shall not need multiply or invert 
such matrices. 

2.2 Calculus on time scales 

Calculus on time scales is a generalization of the standard differential calculus 
and the calculus of finite differences. 

A time scale T is an arbitrary nonempty closed subset of the set R of real 
numbers. In particular T = R, T = hZ for h > and T = q N := {q k , k G N} for 
q > 1 are time scales. We assume that T is a topological space with the relative 
topology induced from R. If to,t\ G T, then [£o,£i]t denotes the intersection of 
the ordinary closed interval with T. Similar notation is used for open, half-open 
or infinite intervals. 

For t G T wc define the forward jump operator a : T— >T by o~(t) := inf {s G T : 
s > <} if t 7^ supT and cr(supT) = supT when supT is finite; the backward jump 
operator p : T— >T by p(t) := sup{s € T : s < t} if t ^ inf T and p(inf T) = inf T 
when infT is finite; the forward graininess function p, : T— »[0,oo) by p{t) := 
o~(t) — t; the backward graininess function v : T— »[0, oo) by v(t) :=t — p(t). 
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If cr(i) > t, then t is called right- scattered, while if p(t) < t, it is called left- 
scattered. If t < supT and <r(t) = t then £ is called right-dense. If £ > inf T and 
p(t) — t, then t is left-dense. 

The time scale T is homogeneous, if and ^ are constant functions. When 
/i = and 1/ = 0, then T = R or T is a closed interval (in particular a half- line). 
When p is constant and greater than 0, then T = /iZ. 

Let T k := {t e T : i is nonmaximal or left-dense}. Thus T fe is got from T 
by removing its maximal point if this point exists and is left-scattered. 

Let / : T^R and t G T fe . The delta derivative of f at t, denoted by / A (i), 
is the real number with the property that given any e there is a neighborhood 
U = (t - 5, t + 5) n T such that 

\(f(a(t)) - f(s)) - f A {t){a{t) - s)\ < s\a(t) s\ 

for all s E U. If f A (t) exists, then we say that / is delta differ entiable at t. 
Moreover, we say that / is delta differ entiable on T fc provided f (t) exists for 
all t G T k . 



Example 2.2. If T = R, then f A (t) = f'{t). If T = KL, then f A (t) 

(qt)-f(t 



f^-fW lfT = g K , then /A (f) _ 



A function / : T— )-R is called rd-continuous provided it is continuous at 
right-dense points in T and its left-sided limits exist (finite) at left-dense points 
in T. If / is continuous, then it is rd-continuous. 

A function F : T— s-K. is called an antiderivative of / : T— >M. provided F A (t) = 
f(t) holds for all t 6 T k . Let a, b G T. Then the delta integral of / on the interval 
[a, 6)t is defined by 

f/(r)Ar:= f f(r)Ar:=F(b)-F(a). 

J a J [a,b)j 

It is more convenient to consider the half-open interval [a, b)j than the closed 
interval [a, b]j in the definition of the integral. If b is a left-dense point, then 
the value of / at b would not affect the integral. On the other hand, if b is left- 
scattered, the value of / at b is not essential for the integral (see Example 12. 3[) . 
This is caused by the fact that we use delta integral, corresponding to the 
forward jump function. 

We have a natural property: 



/ /(r)Ar= f C f(T)Ar+ f /(r)A7 



for any c G (a,b)x- Moreover, if / is rd-continuous, /(<) > for all a < t < b 

b 

and / /(r) At = 0, then / = 0. 
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b b 

Example 2.3. If T = R, then / /(t)At = / /(r)rfr, where the integral on 

a a 

b 

the right is the usual Riemann integral. If T = frZ, h > 0, then J /(r)Ar = 

a 

£-1 

£ f(th)h for a < 6. 

*=t 

2.3 Linear systems on time scale 

Let us consider the system of delta differential equations on a time scale T: 

x A (t) = Ax(t), (1) 
where x(t) £ K n and A is a constant n x n matrix. 

Remark 2.4. If T = M, then ([l) is a system of ordinary differential equations 
x' = Ax. But for T = Z, (|TJ| takes the difference form x(t + 1) — x(t) — Ax(t), 
which can be transformed to the shift form x(t + 1) = (I + A)x{t). Thus to 
compare the definitions and the results stated for delta differential systems in 
the case T = Z with those that were obtained for discrete-time systems in the 
shift form, one has to take this into account. One can easily transform the 
difference form to the shift form and vice versa. 

Proposition 2.5. Equation ([TJ with initial condition x(to) = xo has a unique 
forward solution defined for all t 6 [io,+oo)T- 

The matrix exponential function (at to) for A is defined as the unique forward 
solution of the matrix differential equation X A = AX, with the initial condition 
X(to) = I- Its value at t is denoted by eA(t,to). 

Example 2.6. If T = R, then e A (t,t Q ) = e^'^l If T = hZ, then e A (t,t Q ) = 
(I + A)^)/ h . If T = q N , q > I, then e A (q k t Q ,t Q ) = Utoi 1 + (? ~ 
for k > 1 and to 6 T. 

Proposition 2.7. The following properties hold for every t,s,r £ T such that 
r <s <t: 

i) e A (t,t) = I; 

ii) e A (t,s)e A (s,r) =e A (t,r); 



Let us consider now a nonhomogeneous system 

x A (t)=Ax(t) + f(t) (2) 

where / is rd-continuous. 

Theorem 2.8. Let to G T. System ([2]) for the initial condition x(to) = xq has 
a unique forward solution of the form 

x(t) =e A (t,t )x + f e A (t,a(T))f(T)Ar. (3) 
Jt a 
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3 Positive control systems 



Let ngNbe fixed. From now on we shall assume that the time scale T consists 
of at least n + 1 elements. 

Let us consider a linear control system, denoted by E, and defined on the 
time scale T: 

x A (t) = Ax(t) +Bu{t) (4) 

where t £ T, x(t) £ R™ and u(t) £ R m . 

We assume that the control u is a piecewise continuous function defined on 
some interval [to, £i)t, depending on u, where to S T and ti £ T or tj = oo. We 
shall assume that at each point t £ [to,ti)j, at which m is not continuous, u is 
right-continuous and has a finite left-sided limit if t is left-dense. This allows to 
solve step by step. Moreover, for a finite ti we can always evaluate x{t\). For 
ti being left-scattered we do not need the value of u at t\, and for a left-dense 
ti we just take a limit of x(t) at t\. 

Definition 3.1. We say that system E is positive if for any to £ T, any initial 
condition xo £ R", any control u : [to, *i)t— and any t £ [to,ti]T, the 
solution i of |U satisfies x(t) £ R™ . 

By the separation principle we have the following characterization. 

Proposition 3.2. The system E is positive if and only if e Ait, to) £ R™ x ™ for 
every t, t G T such that t>t , and B £ R" xm . 

The proof is very similar to the proof of the continuous-time case. 

To state criteria of nonnegativity of the exponential matrix, let p, — sup{/x(t) : 
t £ T} and Af := A + I/fl, where //oo means the zero n X n matrix and I/O 
is a diagonal matrix with oo on the diagonal. Thus for T = R, Af is obtained 
from A by replacing the elements on the diagonal by oo, for T = Z, Aj = A + I, 
and for T = q N , A r = A. 

Theorem 3.3. The exponential matrix eA(t, to) is nonnegative for every t, to £ 
T such that t > t if and only if A r £ R+ x ™. 

Proof. "-«=" Assume that A T > 0. If fi(t ) > 0, then A + I/n(t ) > 0. This 
means that e J 4(cr(io), to) = fJ.(to)A + I > 0. If (J,{to) = 0, then for t > to and 
close to to, I + A(t — to) > 0. The last term approximates e Ait, to)- Since 
the exponential matrix is continuous (with respect to t), then also e Ait, to) > 
for t > to and close to to- To achieve nonnegativity of e^(i, to) for all t £ 
T, t > to, we have to use the semigroup property of the exponential matrix: 
e^(t, s)eAis, t) — eA(t, t) for r < s < t and r,s,t £ T. 

Assume that e^t, to) is nonnegative for t, to £ T such that t > to. Suppose 
first that p, > and chose to £ T with jit(to) > 0. Then e^u (to), to) = / + 
fJ>(to)A > 0. This means that also A + I/fJ,(t ) > 0. As it holds for all t £ T 
with n(to) > 0, Af = A + I/p, is nonnegative. If /2 = 0, then T is a standard 
interval. The exponential matrix is then standard e Atyt ~ to \ For t close to to, 
it may be approximated by / + A(t — to). Nonnegativity of the exponential 
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matrix implies that / + A(t — to) > for t > to and close to to- This holds 
only if all elements of A outside the diagonal are nonnegative. Thus again Aj 
is nonnegative. □ 

Corollary 3.4. The system £ is positive if and only if Ay G R + Xn and B G 

jj^n X rn 

Remark 3.5. An nxn matrix with nonnegative elements outside the diagonal 
is called a Metzler matrix. Thus in the continuous-time case, the exponential 
matrix eA(t,to) is nonnegative for every t > to if and only if A is a Metzler 
matrix. In that case the elements on the diagonal may be arbitrary. On the 
other hand, if the time scale T is the set Z of integer numbers, then fi = 1 and 
nonnegativity of the exponential matrix is equivalent to A + / > 0. In that 
case the delta differential equation x A (k) = Ax(k) may be rewritten in the shift 
form as x(k + 1) = (A + I)x(k). Thus the condition A + I > agrees with the 
necessary and sufficient condition of nonnegativity for discrete-time systems of 
the form x{k + 1) = Fx(k), where k G Z (see Q3]). 



4 Controllability 

If £ is a positive system, then for a nonnegative initial condition xq and a 
nonnegative control u, the trajectory x stays in R™. One may be interested in 
properties of the reachable sets of the system. For simplicity we assume that 
the initial condition is xo — 0. Let x(t\, to, 0, u) mean the trajectory of the 
system corresponding to the initial condition x(to) = and the control w, and 
evaluated at time t\. We shall dehne various controllability properties. 

Definition 4.1. Let to,t\ G T, to < t%. The positive reachable set (from 
0) of the system £ on the interval [^Oj^iIt is the set T?.^ '* 1 ' consisting of all 
x(ti,to,0,u), where u is a nonnegative control on \to,ti)T. 
The positive reachable set (from 0) for the initial time to of £ is 

ft+ = U n + M 

ti£T,ii>io 

and the positive reachable set (from 0) of £ is 

The system £ is positively accessible on [£o^i]t if T?.^ '* 1 ' has a nonempty in- 
terior, £ is positively accessible for the initial time to if TZ^ has a nonempty 
interior and £ is positively accessible if 72.+ has a nonempty interior. 
The system £ is positively reachable on [to,ti]j if 7^ 0,tl l — K™^ £ is positively 
reachable for the initial time to if TZ ? = R™ and £ is positively reachable if 
TZ+ = M™ . 
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Remark 4.2. Accessibility was first introduced for nonlinear systems for which it 
is a good substitute of reachability, as the latter is often too restrictive property. 
The same happens for positive systems. Positive accessibility means precisely 
accessibility but with nonnegative controls. 

The following implications follow directly from definitions: 

Proposition 4.3. £ is positively accessible on [to,ti]j £ is positively acces- 
sible for the initial time to £ is positively accessible. 

£ is positively reachable on [toi^ilx S is positively reachable for the initial 
time to £ is positively reachable. 

Positive reachability (on [to,ti]j) implies positive accessibility (on [to,ti]j). 

We have also a useful inclusion: 
Proposition 4.4. If t q < t < t x , then 1Z [t ^ M] C 1Z [ l°' tl] and Tl% C TZ^ . 

Proof. Since we start at xo — 0, to reach x\ € T?.^ '* 1 ' for the initial time To put 
u{t) — for r <E [T ,t ). Then switch to the control that was used to reach X\ 
for the initial time to. □ 

Remark 4.5. Since T may not be homogeneous, in general 7Z+ depends on to- 

Theorem 4.6. Let £ be a positive system and to,ti be elements of T such 
that [to,ti]j consists of at least n + 1 elements. The following conditions are 
equivalent: 

a) £ is positively accessible on [io,^i]T> 

b) £ is positively accessible for the initial time to, 

c) £ is positively accessible, 

d) rank(B, AB, A n ~ 1 B) = n. 

Proof, a) b). This follows from the fact that 7?.+°'*^ is contained in TZ^. 

b) c). This follows from the fact that IZ*! is contained in 1Z+. 

c) => d) . Assume that £ is positively accessible. The reachable set 1Z of £ is the 
set of states that can be reached with controls that are not necessarily nonneg- 
ative. It is clear that IZ^ is contained in 1Z. Moreover, 1Z is a linear subspace 
of 1" (see [2]). Positive accessibility implies that 1Z contains an open subset. 
Therefore 1Z must be equal to W 1 , which means that £ is reachable (controllable) 
from 0. This is characterized by the condition rank(£?, AB, . . . , A n ~ 1 B) — n (see 

d) =>■ a). The condition rank(i?, AB, . . . ,A n 1 B) — n implies reachability of 
£ from on arbitrary interval consisting of at least n + 1 points (see [2]). Ac- 
tually the states can be obtained with the aid of piecewise constant controls 
with at most n — 1 switching at fixed instances. Thus the reachable set IZ' 1 * '* 1 ^ 
on [to,ti]f can be described as the image of a linear map defined on a finite- 
dimensional space. This map restricted to positive controls will give the set 
with the nonempty interior and this means positive accessibility. □ 

To study positive reachability let us introduce a modified Gram matrix re- 
lated to the control system. 
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Definition 4.7. Let M C {1, . . . , m} and i 0) *i S T, to < t\. For each k G M 
let Sk be a subset of [to,ti)j that is a union of finitely many disjoint intervals 
of T of the form [r , ti)t, and let <Sj\/ = {Sk ■ k G M}. By the Gram matrix of 
system corresponding to to, t\, M and Sm we mean the matrix 

W:=W£(M,S M ):=J2 / ^))b k b T k e A {h, a(r)) T Ar. (5) 

Then we have the following characterization: 

Theorem 4.8. Let to,t\ G T. to < t\. System Q is positively reachable on 
[to,ti]i iff there are M C {l,...,m} and the family Sm = {Sk : k G M} of 
subsets of [to,ti]j such that the matrix W = W t * (M,Sm) is monomial. 

Proof. Let x G K" . By e\,. . . , e m we denote the vectors of the standard 

basis in R m . Define control u : [to, >K+ by u(t) = Y^keM u k{ T )^k, where 
Uk{r) = b\"e A (ti, o-(r)) T W -1 x for t G Sk and u k {r) — for t S k - The control 
u is nonnegative and 

e A {h,a(T))Bu(T)AT= ^ / e A {h,a(T))b k u k (T)AT 

= J2 f e J 4(ii,tT(T))6 fe 6je j4 (ii,CT(T)) T VF- 1 5Ar = x. 
fceAf ^ Sfc 

Thus ((!]) is positively reachable on [io,ii]T- 

"=>•" Positive reachability implies that all the vectors ei, . . . , e„ can be reached 
using nonnegative controls. Let us fix some ej. Then there is a piecewise 
continuous nonnegative control u = (u\, . . . , u rn ) on [to,ti)j such that 



E 

3=1 



ti 

eA(ti,a(T))bjUj(T)AT. 



Since all the integrals in the sum are nonnegative, for some fcj the integral 

e j4 (ii,cr(T))6 fei M fei (T)AT 



is an i-monomial vector. Then for every r G [to, ii)T the vector e A (t\, cr(r))6fe i ttfe i (r) 
is either i-monomial or 0. Let £?i be the set all r for which e A (ti, a{T))bk i Uk i { T ) 
is i-monomial. Then for r G Bi the matrix 

(i i , ct(t) )6 fci 6 J. (i i , ct(t)) t 

is diagonal with the only nonzero element at the ith place. The same is true for 
the matrix f B e A (t\, o"(T))6fc 1 6^.e^(ti, <t(t)) t At. This implies that the matrix 

C:=V / e A (i 1 , f r(r))& fe( 6j i e A (t 1 ,cr(r)) T AT 
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is monomial (and diagonal). Let M consist of all fcj for i = 1, . . . ,n. Observe 
that if ki = kj for i ^ j, then Bi n Bj = 0. Define Sk = \J k =k Bi and let 
S M = {S k :kE M}. Then 

C =Y. I e A (t u a(T))b k ble A (t u a(T)) T AT^W^(M,SM), 



fceM 

so W^(M,Sm) is monomial. □ 
Corollary 4.9. If the ordinary Gram matrix 

f ti 

W%= / e A {t u a(T))BB T e A (t 1 ,a(T)) T AT 

J to 

is monomial, then system Q is positively reachable on [to,ti]f. 

Proof. Observe that W% = (M, S M ) for M = {1, . . . , m} and S k = [t , h) r 
for all k e M. Thus positive reachability follows from Theorem 14.81 □ 



Remark 4.10. The condition that is monomial is not necessary for positive 
reachability on [£o,*i]t- Consider the system 



x A = 



-1 1\ 1 1 
1 0/ + U 1 



(6) 



on T = Z. Choose to = and ii = 2. System © is positively reachable on 
[t ,ti] T . Indeed, let M = {1} and Si = [0,2) T . Then 

W =hbj + (I + A)hbJ(I + A) T 



is a monomial matrix. However 

Wh = BB T + (I + A)BB T (I + A) T 



3 3 
3 6 



is not monomial. 



Corollary 4.11. If there exists M C {1, . . . , m} such that the matrix 



is a 



W t ^(M) = J t 1 e A (ti, o~(t))BB t e A (ti, ct(t)) t At is monomial, where B i. 
submatrix of B consisting of column b k ,k € M, then system (j4|) is positively 
reachable on \to,t\\f. 

Proof. Observe that W% (M) = (M, <S M ) where S fe = [t , *i)t for all fc e M. 
Thus positive reachability follows from Theorem 14.81 □ 
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Remark 4.12. The condition that W t 4l (M) is monomial is not necessary for 
positive reachability on [io,*i]T- Let the time scale T = {0} U [1,2] U {3}. 
Consider the system 

The system is positively reachable on [0, 3]t- Indeed, let M = {1} and let 
St = [0,1) T U [2,3) t . Then 

W = [ e A (3,a( T ))BB T e A (3,a(T)) T At 
J[o,i) T 

+ [ e A (3,a( T ))BB T e A (3,a(T)) T AT 

J[2,3) T 

o o Wi oWi o 

e- 2 J + [O Oj ~ {0 e~ 2 

is monomial. Observe that we remove here the points t with fi(t) = 0. This is 
essential in order to get a monomial matrix. To calculate the full Gram matrix 
we have to add to W the following matrix 



/ e A (3,a(T))BB T e A (3,a(T)) T dT. 



Its off-diagonal elements are equal to f?(3 — r)e 2 ^ 3 T ^dr. Since they are posi- 
tive, W t 1 (M) is not monomial. 

For T = R we obtain very restrictive conditions for positive reachability. 

Proposition 4.13 {[6 ). Let T = R and i ,*i £ R, i < *i. System @) 
is positively reachable on [to,ti] iff A is diagonal and B contains an n x n 
monomial submatrix (so m > n). 

Proof. "<J=" Let B denote the monomial submatrix of B and let the indices of 
columns of B form the set M. Then BB T is a diagonal matrix with all the 
diagonal elements being positive and so is 

/•ti 

W&(M) = / e A (t 1 ,a(T))BB T e A (t 1 ,a(T)) T AT. 

Thus W^(M) is monomial, so system Q is positively reachable by Corol- 
lary [4jT] Observe that the proof of this implication works for all time scales. 

"=>" Assume that the system is positively reachable on [to, t{\. From Theo- 
rem 14.81 it follows that for some set M and some family Sm the Gram matrix 
W = W\l (M, Sm) is monomial. Let jth column of W be i-monomial. Then for 
some k £ M and for r from some subinterval of [io>*i) the jth column of the 
matrix e A ( tl ~ T " l bkb'£(e A ( tl ~ T ' > ) T is i-monomial. Let c(r) = e A ( tl ~ T ' ) bk- Since the 
jth column of the matrix c(t)c(t) t is z-monomial, then c(r) must be j-monomial 
and eventually i — j. This means that at least one column of e A ^ ll ~ T ^ must 
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be i-monomial. As the exponential matrix is invertible such a column must be 
unique. This implies that bk is monomial. Moreover the i-monomial column of 
e A(ti-r) mus t De j^g Hh column. Otherwise we would get on the diagonal of 
the exponential matrix for all r from some interval, which is impossible. Thus 
e A(ti-T) j g ^^0^1 on SO me interval, which means that A is also diagonal. Now 
to get all n monomial columns in W we need n different monomial column bk- 
Thus B contains an n x n monomial submatrix. □ 

Remark 4.14. The statement of Proposition l4.13l holds also for T = [a, b] and for 
T being a closed half-line. However it does not hold for disjoint union of closed 
intervals. The example given in Remark 14.121 may be considered on a bigger 
time scale: T = [—1, 0] U [1, 2] U [3, 4]. Neither A nor B satisfy the requirement 
given in Proposition ^. 131 However the system is positively reachable on [0, 3]t- 

Corollary 4.15. For T = M., system (j4|) is positively reachable on some [to,ti] 
iff (j4]) is positively reachable on any interval [to,ti]. 

Corollary 4.16. For T = R, system (U) is positively reachable on some [to,ti] 
iff (j4]) is positively reachable. 

For discrete homogeneous time scales the conditions for positive reachability 
are much less restrictive. 

Proposition 4.17. Let T = fiZ for a constant /i > 0. Let to e T and t\ = 

to + k[i for some k G N. System ([4J is positively reachable on \po,ti]j iff the 
matrix [B, (L + [iA)B, ...,(/ + /iA) k ~ 1 B] contains a monomial submatrix. 

Proof. Observe that x(h) = E^oEjLil 7 + nA)%Uj{k - 1 - i). If 

(I + /j,A) l bj = je s for some 7 > 0, then setting Uj(k — 1 — i) = 1/7 and all other 
components and values at different times putting to we get x(t\) = e s . This 
means positive reachability on [fo,ii]T- 

"=>" By Theorem 14.81 positive reachability implies existence of a set M and 
subsets Sk of [to, ^i] for k 6 M such that the matrix 

W=J2 [ eA(ii,a(T))6 fe 6^(ii,a(T)) T Ar 

is monomial. Moreover 

/ e A (t u o(T))b k ble A (t u a(T)) T AT = 
Js k 

^2(1 + MA) (tl -'-^& fe &£((J + nA^-t-rtl^fn. 
test 

This implies that for every i = 1, . . . , n there are k £ M , t E S k and < j < n 
such that the jth column of (/ + ^A)^ 1- * - ^ ^ l b k bl((I + ^A)^ 1 -*-^/^ is 
i-monomial. This means that the column (/ + iJ,Ap tl ~ t ~ li '' IJ, bk is i-monomial. 
But this column is one of the columns of the matrix [B, (I + fiA)B, ...,(/ + 
liAf^B]. □ 
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If k > n then it is enough to consider powers of I + [iA only up to n- 1. 



Proposition 4.18. Let T = /iZ, k > n, to G T and t-y = to + kfi. System (U|) is 
positively reachable on [to, ti]j iff the matrix [B, (I + [j,A)B, ...,(/ + iiA) n ~ l B\ 
contains a monomial submatrix. 

For /i = 1 this was shown in ]S]. Then Q may be rewritten in a more 
familiar form x(t + 1) = (I + A)x(t) + Bu(t). The proof for ^ / 1 is very 
similar. 

Proposition 14. 17l may be extended to nonhomogeneous discrete time scales. 

Proposition 4.19. Assume that fi(t) > for all ( e T, ^ £ T and t\ = a k (to). 
System ([4]) is positively reachable on [to,ti]f iff the matrix 

[B, (I + n(<r{to)A)B,(I + ^(a 2 (t )A)(I + n(a(t )A)B, 

(I + fi^- 1 (t ) A) ...(/+ fi(a(t )A)B] 

contains a monomial submatrix. 

The proof is similar to the proof of Proposition 14.171 but we have to take 
into account that the exponential matrix is no longer a power of I + [lA for a 
constant p, but rather a product of such terms with possibly different values of 
[i. This criterion may be used for systems on T = q N . 

Remark 4.20. Proposition ^. 18l cannot be extended to discrete nonhomogeneous 
time scales. Let us consider the following example: T = {0, 1, 2, 4}, n = 2, 

*-(?4)--G 

Let to — and t\ = 4, so k = 3. The matrix 

[B,(I + f x(l)A)B,(I + f,(2)A)(I + ^l)A)B] = (l \ \ 

evidently contains a monomial 2x2 submatrix. But this is not true for the 
matrix [B, (I + [i(\)A)B\. Thus to reach monomial vectors we may need more 
than n jumps. 
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